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ABSTRACT 

A probabilistic approach is described for aeropropulsion system as- 
sessment. To demonstrate this approach, the technical performance of a 
wave rotor-enhanced gas turbine engine (i.e. engine net thrust, specific 
fuel consumption, and engine weight) is assessed. The assessment 
accounts for the uncertainties in component efficiencies/flows and 
mechanical design variables, using probability distributions. The results 
are presented in the form of cumulative distribution functions (CDFs) 
and sensitivity analyses, and are compared with those from the tradi- 
tional deterministic approach. The comparison shows that the probabi- 
listic approach provides a more realistic and systematic way to assess an 
aeropropulsion system. 

INTRODUCTION 

The need to provide cost-effective aeropropulsion technology is 
critical in today’s environment. Since the beginning of the 1990s, civil 
and military organizations have been attempting to reduce costs like never 
before. Demands are strong for aircraft engines with reduced life-cycle 
costs, emphasizing durability for longer intervals between overhaul, high 
reliability for the lowest possible unscheduled maintenance, and improved 
maintainability for fast, simple maintenance. In response to such demands 
of their highly competitive marketplace, the engine manufacturers must 
focus on product reliability, maintainability, and most definitely 
affordability. At the same time they must reduce the product develop- 
ment cycle due to shrinking development budgets and resources. Rapid 
turn around time while investigating new design concepts or technolo- 
gies is critical to be competitive within the aerospace industry. At the 
Propulsion System Analysis Office (PSAO) of NASA Glenn Research 
Center (GRC), our mission is to assess the performance potential and 
economic benefits of advanced and unconventional propulsion systems 
for a broad spectrum of subsonic, supersonic, and hypersonic aeronauti- 
cal vehicles. These assessments provide the basis for NASA’s future 
aeropropulsion program directions and technology investment decisions. 
As the lead NASA Center in aeropropulsion, GRC’s primary mission is 
aeropropulsion research and technology. The technology is transferred 
to the aviation industry to help maintain U.S. leadership in the world s 
aviation market. 

In view of the challenges facing the engine manufacturers today, it 
is obvious that critical decisions must be made in the early stages of 
engine development where available design freedom can best achieve 
better performance and cost is relatively low, as shown in Fig. 1 . As such, 
the role of system assessment, performed in the early stages of an engine 


development program, becomes very critical to the successful develop- 
ment of new aeropropulsion systems. A reliable system assessment not 
only helps to identify the best propulsion system concept amongst sev- 
eral candidates, it can also identify high payoff technologies worthy of 
pursuit to decision-makers. This is particularly important for advanced 
aeropropulsion technology development programs which require enor- 
mous amount of resources, such as the Pulse Detonation Engine (PDE) 
and Supersonic Through-Flow Fan (SSTF) concepts being investigated 
at GRC. 

In the current practice of deterministic or point-design approach, 
uncertainties of design variables are either unaccounted for or accounted 
for by the safety factors. This could often result in an assessment with 
unknown and unquantifiable reliability. Consequently, it fails to pro- 
vide additional insight into the risks associated with new technologies, 
which are often needed by the decision-makers to determine the benefit 
and return-on-investment of a new aircraft engine. In this paper, an al- 
ternative approach based on probabilistic method is described for a com- 
prehensive assessment of aeropropulsion system. The statistical 
approach quantifies the design uncertainties inherent in a new 
aeropropulsion system and their influences on engine performance. It 
provides an analytical framework that allows an engine developer to 
improve engine performance by determining the necessary design mar- 
gin, the parameters impacting the uncertainty in performance, and ways 
to reduce the impact of uncertainty. As such, it enhances the reliability 
of a system assessment. An assessment of a wave rotor-enhanced tur- 
bofan engine is performed to demonstrate the methodology. The as- 
sessment accounts for the uncertainties which occur in component 
efficiencies/flows, and mechanical design variables, using probability 
distributions. 

OBJECTIVE 

The objective of the current work is to demonstrate the application 
of probabilistic approach and its feasibility for aeropropulsion system 
assessment. 

NUMERICAL EXAMPLE FOR DEMONSTRATION 

A wave rotor-enhanced turbofan engine is chosen to demonstrate 
the probabilistic approach. Wave rotor is a device that utilizes unsteady 
wave to compress air in a single device. Enhancement of a gas turbine 
engine with a wave rotor can improve specific power and reduce specific 
fuel consumption. This engine was analyzed previously using the tradi- 
tional deterministic approach, and the results are reported (Jones and 
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Figure 1.— Design process paradigm (Aerospace System Design Laboratory, 1996). 


Welch, 1996). The intent of the current work is not to challenge the 
previous results. Rather, it is to demonstrate the application and feasibil- 
ity of probabilistic approach for aeropropulsion system assessment. 

Affordability and retum-on-investment of a new aeropropulsion 
system are strongly influenced by cost, specific fuel consumption, 
engine thrust and weight. This demonstration focuses on the technical 
aspect of engine performance, i.e., specific fuel consumption, engine 
thrust and weight. 

The current results are compared with those calculated previously 
using the deterministic approach (Jones and Welch, 1996). In addition, 
the baseline engine used for the comparison with the wave rotor-enhanced 
engine (Jones and Welch, 1 996) is also used for the current comparison. 

ANALYSIS APPROACH AND PROCEDURES 

The approach taken in this effort is to combine thermodynamic cycle 
analysis embedded in the computer code NEPP (NASA Engine Perlor- 
mance Program, Klann and Snyder, 1994), engine weight analysis 
embedded in the computer code WATE (Onat and Klees, 1 979), and fast 
probability integrator (FP1, Southwest Research Institute, 1995). A sche- 
matic of the integrated approach is shown in Fig. 2. 

The computer code NEPP is used to calculate engine net thrust and 
specific fuel consumption. It is a one-dimensional steady state thermo- 
dynamic cycle analysis code which allows the user to model virtually 
any kind of gas turbine engine cycle through the use of components which 
can be placed in any order to create the desired cycle. The engine weight 
is estimated by the WATE code. These analyses are performed at the sea- 
level static condition. The role of FPI is to perform probabilistic analysis 
utilizing the results generated by NEPP and WATE. The results are 


presented in the form of cumulative distribution functions (CDFs). In 
addition, FPI is used to perform sensitivity analysis to rank design vari- 
ables in order of their impact on a specific response variable. The FPI 
code was developed under contract with NASA Glenn Research Center. 

Sensitivity values could be ± in nature. A positive value indicates 
that the response variable will decrease as the design variable increases 
and vice versa. Variable with the highest absolute sensitivity value is 
defined to be the most influential variable. Variable with the second highest 
absolute sensitivity value is second most influential variable and so on. 
This ranks the design variables in the order of their influences on the 
response variable. The sensitivity information thus obtained from FPI is 
very useful from the design point of view. For example, reliability in 
design can be improved when uncertainties in the most influential vari- 
ables are reduced. Those design variables that do not have significant 
influences deterministically could nevertheless have strong influences 
on the response scatter if these design variables have large uncertainties. 
Weak design variable with large uncertainties may have probabilistic 
sensitivity factors more important than strong design variables with small 
uncertainties. Unlike deterministic analysis, sensitivity factors in proba- 
bilistic analysis arc functions of both the deterministic sensitivity and 
the uncertainty (characterized by the standard deviation). 

The procedures for the probabilistic analysis are as follows: 

( 1 ) Identify design variables with uncertainties (i.e. identify the risk 
elements). 

(2) Quantify the uncertainties with probability distributions, based 
on expert opinion elicitation, historical data, or benchmark/ 
prototype testing, etc. 
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Figure 2.— Fast probability integration input/output schematic. 


(3) Perform perturbation for the selected set of values (mean and 
standard deviation) of the design variables to generate response 
variables. For this work, NEPP and WATE are used to perturb 
specific fuel consumption (SFC), engine thrust and weight. 

(4) Perform probabilistic analysis to compute CDFs of the response 
variables. For this work, FPI is used to compute the CDFs and 
the corresponding sensitivities of the response variables and rank 
the design variables in the order of their influences. 

PROBABILISTIC SIMULATION 

There are a number of approaches available for obtaining probabi- 
listic response for a given set of independent primitive variables. One 
approach that is commonly used is called Monte-Carlo simulation tech- 
nique. While this approach provides almost exact solutions, it is fairly 
expensive and time-consuming approach. In this technique, randomly 
selected values of the input variables, which are based on their known 
probabilistic distributions, are used to deterministically compute the value 
of the response variable. This has to be repeated usually several hun- 
dreds or even thousands of times to build the response probabilistic char- 
acteristics. In essence, this technique requires a large number of 
simulations to generate CDF’s of output variables. Although, inherently 
simple, the large number of computer runs required to obtain a reason- 
ably accurate CDF of output variables becomes its obvious disadvan- 
tage. NASA-Glenn Research Center has been involved in developing 
efficient probabilistic methods for more than a decade. As a result of this 
research initiative, fast probability integration (FPI) algorithms were 
developed (Aerospace System Design Laboratory, 1 996) to solve a large 
class of engineering problems. 

Let us say that there are n random design variables in a problem and 
that we want to use probabilistic analysis to compute the probability of 
occurrence of a certain response function 

Z(X) = Z(X { ,X 2 ...X n ) ( 1 ) 

where Z represents the response variable and X represents the random 
variable. Our aim is to compute the probability that Z will have a value 


less than or equal to a given magnitude Z 0 . To achieve this goal, the 
performance function, which describes how the mechanics of the system 
behaves, can be cast as a limit state tunction g(30, which can be 
described as 

g(X) = Z(X)-Z 0 ( 2 ) 

Traditionally, the limit state function g has been defined in such a 
way that g = 0 represents a boundary where g < 0 represents tailed region 
and g > 0 represents safe region. Here the objective would be to compute 
P[g(X) < 0]. Generally speaking, Z Q describes a limit indicating failure; 
g(X) is called a failure function. 

Given the joint probability density function f x (x) of the limit state 
function g(x ), we can formulate the limit state probability P[g < 0] as 

P f = F{g{X)<0\ = y..\f x (x)dx ( 3 ) 

a 

where a describes the domain of integration. This multiple integration 
is, in general, very difficult to integrate analytically. However, FPI has 
been found to be an excellent tool to evaluate Eq. (3) efficiently and 
accurately. 

DESIGN AND RESPONSE VARIABLES FOR THE 
NUMERICAL EXAMPLE 

The design variables are: 

fan efficiency 

low pressure compressor (LPC) efficiency 
high pressure compressor (HPC) efficiency 
wave rotor pressure ratio 
high pressure turbine (HPT) efficiency 
high pressure turbine (HPT) inlet temperature 
low pressure turbine (LPT) efficiency 
bleed How 

turbine disk material strength 
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TABLE 1. — Wave Rotor-Enhanced Turbofan Engine 
Design Variables wilh Uncertainties 


Design variable 

Deterministic 

Probabilistic Approach 


approach 

(Jones and Welch, 1996) 

Mean 


■ 


0.91 

0.91 

±0.02 

Normal 

LPC efficiency 

0.88 

0.87 

±0.02 

Nonna! 

HPC efficiency 

0.85 

0.87 

±0.02 

Normal 

Wave rotor pressure ratio 

1.15 

1.13 

±0.02 

Normal 

HPT efficiency 

0.89 

0.88 

±0.02 

Normal 

HPT inlet temperature 

3200 R (1778 K) 

3200 R (1778 K) 

±100 R (±56 K) 

Nonnal 

LPT efficiency 

0.93 

0.91 

±0.02 

Normal 

Bleed flow, percent 

19.5 

19.0 

±1.0 

Normal 

Turbine disk material strength 

100 ksi 
(690 MPa) 

100 ksi 
( 690 Mpa) 

±10 ksi 
(69 Mpa) 

Weibull 


Other Design Variables 

Design variable 

Deterministic approach 

Probabilistic approach 

Inlet flow 

2800 lb/s (1270 kg/s) 

Inlet recovery 

1.00 

Inlet temperature 

545.7 R (303 K) 

Bypass ratio 

7.00 

Fan pressure ratio 

1.59 

Fan corrected flow 

2875 lb/s (1 304 kg/s) 

LPC pressure ratio 

1.55 

HPC pressure ratio 

15.8 

Wave rotor temp, ratio 

1.9! 


TABLE 2. — Performance of the Wave Rotor-Enhanced Turbofan Engine 
Comparison of Results between Probabilistic and Deterministic Approaches 


Engine performance 

Deterministic approach 
(Jones and Welch, 1996) 

Probabilistic approach 

Specific fuel consumption 

0.304 lb/hr/lb 
(0.0310 kg/hr/N) 

< 0.304 Ib/hr/lb (0.03 10 kg/hr/N) 1 0% probability 

< 0.309 Ib/hr/lb (0.0315 kg/hr/N) 50% probability 

< 0.3 16 lb/hr/lb (0.0322 kg/hr/N) 95% probability 
Median value = 0.309 Ib/hr/lb (0.0315 kg/hr/N) 

Engine net thrust 

89470 lbs 
(398 kN) 

> 89470 lbs (398 kN) 60% probability 

> 901 14 lbs (401 kN) 50% probability 

> 85997 lbs (383 kN) 95% probability 
Median value = 901 14 lbs (401 kN) 

Engine weight 

‘2 11 20 lbs 
(9580 kg) 

<211 20 lbs (9580 kg) 66% probability 
< 20952 lbs (9504 kg) 50% probability 
<21650 lbs (9820 kg) 95% probability 
Median value = 20952 lbs (9504 kg) 


Note: *AU engine weight calculations include the weight of the wave rotor, estimated to be 1 650 lbs 
(748 kg) (Jones and Welch, 1996), 


The assumed mean values and standard deviations of these design 
variables are shown in Table 1 . These variables are assumed to be 
independent and have normal distributions, with the exception of the 
turbine disk material strength which has Weibull distribution. For many 
engineering problems, design variables with small variability have gen- 
erally been seen to have normal distributions. Material strength is gener- 
ally characterized with Weibull distributions. With upper and lower limits 
for these design variables based on ±2 standard deviations around the 
mean value for each design variable, 95.4 percent of the expected uncer- 
tainty is captured. 


The response variables are: 

engine net thrust 

specific fuel consumption (SFC) 

engine weight 

RESULTS AND DISCUSSION 

As mentioned earlier, it is critical to assess the reliability of a new 
aeropropulsion system because of inherent design uncertainties. The 
current demonstration focuses on the technical aspect of engine perfor- 
mance, i.e., specific fuel consumption, engine thrust and weight. The 
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Note: + SFC decreases as design variable increases. 

— SFC increases as design variable increases. 

Figure 4.— Sensitivity of specific fuel consumption. 


results are presented in the form of cumulative distribution functions 
(CDFs) and sensitivity analyses. A CDF gives a relationship between a 
value up to certain magnitude of a response variable and the probability 
of its occurrence. 

The results are summarized and compared with those based on the 
deterministic approach, in Table 2. 

In reference to the SFC obtained previously using the deterministic 
approach, 0.304 lb/hr/lb (0.03 10 kg/hr/N), the current result shows that 
the cumulative probability for that to occur is only 10 percent. In 
other words, the confidence level to achieve a SFC of 0.304 Ib/hr/lb 


(0.0310 kg/hr/N) or better is only 10 percent. The probability is much 
higher to obtain a SFC of 0.3 1 6 Ib/hr/lb (0.0322 kg/hr/N) or better, about 
95 percent. This is shown in Fig. 3. 

The sensitivity of SFC to the nine design variables is shown in 
Fig. 4. It shows that the most influential design variable for SFC is high- 
pressure turbine (HPT) inlet temperature. It implies that to improve the 
SFC of this engine, the biggest payoff is to decrease the HPT inlet tem- 
perature and control its scatter (uncertainties). By reducing the mean 
temperature from 3200 R ( 1 778 K) to 3 1 60 R ( 1 755 K), the SFC reduces 
from 0.309 lb/hr/ib (0.0315 kg/hr/N) to 0.306 Ib/hr/lb (0.031 2 kg/hr/N) 
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Figure 5. — PDF of wave rotor-enhanced turbofan engine SFC. 



at 50 percent probability level. This result is also shown in Fig. 3. In 
addition, by reducing the uncertainty of the HPT inlet temperature from 
± 1 00 R (56 K) to ±50R (28 K), the scatter range of the SFC is reduced by 
30 percent. This is indicated by the probability density function (PDF) 
shown in Fig. 5. The influences of HPT and LPT efficiencies are moder- 
ate. Other design variables such as LPC and HPC efficiencies, and wave 
rotor pressure ratio have minimal influences. As expected, the turbine 
disk material strength has no influence on SFC. 

For the engine net thrust, the result shows that the probability to ob- 
tain an engine net thrust of 89470 lbs (398 kN) or higher, which was re- 
ported (Jones and Welch, 1996) based on the deterministic approach, is 
fairly good — about 60 percent (1.0 to 0.4). This is shown in Fig. 6. Also, 
even with the current assumption of design uncertainties, it is very likely 


that the current engine will outperform the baseline engine in net thrust — 
about 90 percent probability. This is also shown in Fig. 6. The sensitivity 
results for engine net thrust is shown in Fig. 7. Again, the result shows that 
the most influential design variable for engine net thrust is HPT inlet tem- 
perature, and the turbine disk material strength has no impact. 

For the engine weight, the likelihood of obtaining an engine weight 
of 21 120 lbs (9580 kg) or less, calculated deterministically (Jones and 
Welch, 1996), is fairly good — about 66 percent. However, with the cur- 
rent assumption of design uncertainties it is very unlikely that the cur- 
rent engine will outperform the baseline engine in engine weight — only 
about 10 percent probability. These results arc shown in Fig. 8. The sen- 
sitivity results for engine weight is shown in Fig. 9. As expected, the 
result shows that the turbine disk material strength has significant 
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Figure 7. — Sensitivity of engine net thrust. 


strength 



impact on the engine weight. The material strength, together with HPT 
inlet temperature and HPC efficiency, are the three most influential design 
variables. The rest of the design variables have moderate influences. 

Overall, the results show that component uncertainty can have a 
significant impact on engine performance. For the current example, the 
HPT inlet temperature dominates the influence on SFC and engine thrust. 
Its impact on engine weight is more moderate, since other design variables 
such as turbine-disk material strength and HPC efficiency also have 
significant impact on engine weight. Decreasing the HPT inlet tempera- 
ture will have favorable impact on SFC but unfavorable impact on 
engine thrust and weight, and vice versa. If SFC receives too much 


emphasis, engine thrust and weight suffer. However, increasing the 
turbine-disk material strength and/or component efficiencies can com- 
pensate the unfavorable impact on engine weight. Discovering these 
relationships is not peculiar to using probabilistic approach. Rather, proba- 
bilistic approach helps engine developers visualize and make trades of 
design margins. 

More over, the results show that the current integrated probabilistic 
approach (NEPP + WATE + FPI) not only calculates the SFC, engine 
thrust and weight, but also determines the probability of their occur- 
rences. As such, the probabilistic assessment provides additional insight 
into the risks associated with new technologies, which makes it easier 
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for the decision-makers to determine the benefit and return-on- 
investment of a new aeropropulsion system. In addition, the approach 
ranks the relative importance of the design variables as to their influ- 
ences on the engine performance, in the form of sensitivity factors. High 
sensitivity factor indicates an area to be focused for improving the 
engine performance and its reliability. 

The current work addresses the application of probabilistic approach 
to assess SFC, engine thrust and weight. Similarly, the approach can also 
be used to assess the other aspects of aeropropulsion system performance, 
such as cost, acoustic noise, and emissions, etc. 

CONCLUSIONS 

Based on the comparison of the results between the probabilistic 
and deterministic approaches, the following conclusions are made: 

(1) The probabilistic approach provides a more realistic and sys- 
tematic way to assess an aeropropulsion system, because it 
accounts for uncertainties in the design variables. 

(2) The results from probabilistic assessment are more credible and 
reliable, because it incorporates the ‘past lessons learned' (i.e., 
expert opinions, historical data, etc.) to quantify the risks. In 
addition, the likelihood of repeating past mistakes will be 
minimized. 

(3) The probabilistic approach allows the decision-makers to detect 
problems early before they become critical. As such, risks can 
be mitigated accordingly and resources (time, funding, etc.) can 
be used more wisely. 

(4) Probabilistic assessments provide decision makers with a tool 
that allows them to assign priorities to needed technological 


developments and thus increase the likelihood that R&D invest- 
ments will have high payoffs. 

(5) The probabilistic assessment provides additional insight into the 
risks associated with new technologies, which makes it easier 
for the decision-makers to determine the benefit and retum-on- 
investment of a new aeropropulsion system. 
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